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George C. Soulas 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Abstract 

The equations for rigorously calculating the particle flux and density surrounding a cylindrical aperture in the 
free molecular flow regime are developed and presented. The fundamental equations for particle flux and density 
from a reservoir and a diffusely reflecting surface will initially be developed. Assumptions will include a Maxwell- 
Boltzmann speed distribution, equal particle and wall temperatures, and a linear flux distribution along the 
cylindrical aperture walls. With this information, the equations for axial flux and density surrounding a cylindrical 
aperture will be developed. The cylindrical aperture will be divided into multiple volumes and regions to rigorously 
determine the surrounding axial flux and density, and appropriate limits of integration will be determined. The 
results of these equations will then be evaluated. The linear wall flux distribution assumption will be assessed. The 
axial flux and density surrounding a cylindrical aperture with a thickness-to-radius ratio of 1 .25 will be presented. 
Finally, the equations determined in this study will be verified using multiple methods. 

Nomenclature 

a y-intercept in non-dimensional wall flux equation 

A area 

b slope in non-dimensional wall flux equation 

c mean thermal speed of particles within the reservoir 

f(v) Maxwell-Boltzmann space speed distribution function 

n particle density 

L distance between the surface normals of dA 0 and dA c 

N particle flow rate 

P particle transmission probability through a cylindrical aperture 

R aperture radius 

R p projected circle radius 

r radial location 

r 0lI radial limit of integration for the reservoir’s contribution to the downstream flux and density 

r oln radial limit of integration for the reservoir’s contribution to the downstream flux and density 

t aperture plate thickness 

T far-field angular distribution function for particles downstream of a cylindrical aperture 

v particle speed 

z axial location 

z wII axial limit of integration for the wall’s contribution to the downstream flux and density 

z wIIn axial limit of integration for the wall’s contribution to the upstream flux and density 

Greek Symbols: 

a angle from dA 0 surface normal to L 

P angle 
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y angle 

F particle flux 

F 0 reservoir or wall flux per steradian 

5 angle from dA t . surface normal to L 

0 polar angle in a spherical coordinate system 

k far-field angular distribution function variable 

9 azimuthal angle 

9 omax angular limit of integration for the reservoir’s contribution to the downstream flux and density 

(Poii angular limit of integration for the reservoir’s contribution to the downstream flux and density 

(p w ii angular limit of integration for the wall’s contribution to the downstream flux and density 

(p w ii n angular limit of integration for the wall’s contribution to the upstream flux and density 

(p wt angular limit of integration for the wall’s contribution to the flux and density 

co solid angle 


Subscripts and Superscripts: 

non-dimensionalized value (dimensions are divided by R, densities by n 0 , and fluxes by r o ) 
c collection site 

f final location for the limit of integration 

far far-field variable 

i initial location for the limit of integration 

o reservoir 

up upstream reservoir 

w aperture wall 


I. Introduction 

Accurately predicting the density and flux of particles emanating from a circular aperture of finite thickness 
(i.e., a cylindrical aperture) is of significant interest for ion thrusters. Ion thrusters create ions within a discharge 
chamber and accelerate these ions through perforated grids composed of thousands of cylindrical apertures to 
produce thrust (Refs. 1 and 2). Xenon is the propellant used for present spacecraft ion thruster applications. One 
consequence of beam ion generation is that unionized gas can escape from the thruster discharge chamber through 
the cylindrical apertures of the grids. Predicting the downstream unionized gas density distribution is important for 
determining its impact on the thruster and surrounding spacecraft environment. For example, slow-moving ions are 
created with beam ions and the escaping unionized gas due to a process known as charge-exchange. These slow- 
moving ions can either impact the thruster’s accelerator grid, causing erosion and limiting thruster service life, or 
they can impact spacecraft surfaces. 

Ion thrusters operate at neutral densities that are so low that particles travel from surface to surface without 
colliding. For example, the neutral-to-neutral mean free paths are greater than 34 cm within a state-of-the-art 
discharge chamber, which is greater than the radius of the chamber (Ref. 3). This rarified gas flow regime is often 
referred to as free molecular flow. Solving for the downstream particle density distribution requires an 
understanding of how the thickness of a cylindrical aperture can affect this downstream distribution. 

Many ion thrusters utilize two grids to accelerate ions, and these include an upstream screen grid and a 
downstream accelerator grid. So, for a unionized gas to escape the discharge chamber, it must travel through two 
apertures, which presents a difficult problem to solve. Fortunately, though, the upstream screen grid aperture is 
typically much larger and thinner than the accelerator grid aperture. As a result, the angular distribution of particles 
exiting the accelerator aperture is well modeled by assuming a single accelerator aperture (Ref. 3). 

Many past studies have considered this problem of the free molecular flow through a cylindrical aperture. 
Clausing was the first to determine the downstream angular distribution pattern of particles emanating from a thin 
cylindrical tube under free molecular flow conditions (Refs. 4 and 5). Flowever, his distribution equation was only 
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appropriate for far-field calculations of particle fluxes and densities, as explained by Reynolds and Richley (Ref. 6). 
Richley and Reynolds determined the flux of particles within converging and diverging cylindrical tubes (Refs. 6 
and 7). And Robertson conducted similar studies to determine particle densities as well (Ref. 8). The results of both 
investigations can be extended axially downstream of the cylindrical aperture to determine near- and far-field 
densities and fluxes. Unfortunately, the resulting equations are not accurate at radial locations greater than the radius 
of the aperture, which imposes a severe limitation on their use. Another method of determining particle density and 
flux distributions is to use a Monte Carlo simulation technique (Ref. 9). Such a technique allows for the quick 
analysis of complex geometries (see, for example. Ref. 10). However, accurately analyzing downstream regions 
where particle densities and fluxes have decreased significantly is computationally difficult. In addition, this 
technique does not lend itself well to ion thruster grids which have thousands of cylindrical apertures. 

In this study, equations will be developed for rigorously calculating the flux and density of particles 
surrounding a cylindrical aperture in the free molecular flow regime. The fundamental equations and their 
assumptions necessary to determine particle densities and fluxes will be presented first. These equations will then be 
applied to a cylindrical aperture to determine the particle flux distribution from the cylindrical aperture walls. With 
this information, the equations for particle density and flux surrounding a cylindrical aperture will be developed. 
Finally, the resulting equations will be applied to a typical ion thruster aperture and compared to that of Clausing in 
the far-field (Ref. 4). 


II. Fundamental Equations for Particle Flux and Density 

The particle flux and density distributions surrounding a cylindrical aperture were determined using the 
following assumptions: 

1 . Particles upstream of a cylindrical aperture reside within a reservoir and exhaust into a vacuum through the 
cylindrical aperture; 

2. Equilibrium conditions exist, so that the rate of particles leaving the reservoir is equal to the rate of 
particles fed into the reservoir; 

3 . The particle density is sufficiently low that free molecular flow conditions exist within the aperture; 

4. All wall reflections are diffuse; 

5. Particle adsorption onto the walls and surface diffusion are negligible; 

6. Wall temperatures are equal to the particle temperature; and 

7. Particles within the reservoir have a Maxwell-Boltzmann speed distribution. 

The geometric relation between a differential element from a particle emission site (i.e., a wall or the upstream 
reservoir) and a differential element from a downstream collection site is shown in Figure 1. 


Collection Site 


dA, T, n 



dA 0 , r 0 , n 0 
Emission Site 


Figure 1. — Geometric relation 
between a particle emission 
site and the downstream 
collection site. 
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The variable a is the angle between the surface normal of the emission site to the collecting site’s differential surface 
area, dA, and 5 is the angle between the surface normal of dA to the differential emitting area, dA 0 . The variable L is 
the distance between the emitting and collecting differential areas. The particle flow rate, N , from a differential 
emission area to a downstream differential collection area is given by: 

d 2 N = f 0 (Aj-dA 0 •cos(a)-dot (1) 


Here, f 0 is the emitted particle flux per steradian from differential area dA 0 and is a function of A 0 . 
The downstream solid angle, dco, is given by: 


doj = 


_ dA • cos(5) 


(2) 


By substituting Equation (2) into Equation (1) and dividing both sides by dA, the following equation is obtained: 


^ dN^ 


v dAy 


■ dr.r 0 (A 0 ). cos (“>r s < 5 > .dA 0 


( 3 ) 


Here, T is the flux to the collection site. The emitted particle flux, T 0 , can be evaluated as a function of the emitted 
particle flux per steradian by noting that for diffuse surface emission, r o is independent of emission angle, so that: 


2-71 2 

T 0 = | Jf D •cos(9)-sin(9)-d0-d(|) = Jt-f 0 (4) 

o o 

Note that particles from the upstream reservoir and reflected from a wall can both be assumed to be leave 
differential area dA 0 diffusely. Substituting this result into Equation (3) and integrating both sides yields our 
fundamental equation for particle fluxes: 

r A.Jr 0 (Aj.f2!W^®.dA 0 (5 ) 

* L L ' 

Note that if the emission site is the reservoir, T 0 is independent of A 0 and equal to: 



(6) 


where c is the mean speed of the particles within the reservoir and n 0 is the reservoir particle density. 

The downstream particle density, n, from a reservoir of particles can be solved by starting with the Maxwell- 
Boltzmann phase space speed distribution function (Ref. 11). Because equilibrium conditions were assumed, the 
speed distribution function, f(v), is independent of position (i.e., particle velocity is isotropic). So, the particle 
density downstream of a reservoir is: 


n = n ° • JJJ f ( v ) d3v = n ° • J jf(v)-dco 0 -dv = n 0 jf(v)-dv- J 


dff> n 


0 co„ 


Here, n 0 is the particle density within the reservoir and the differential solid angle, dco 0 , is given by: 


dro n = 


dA 0 -cos(a) 


( 7 ) 


( 8 ) 
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The first integral on the right-hand side of Equation (7) can be shown to be 1/(4 -ti), so substituting this and 
Equation (8) into Equation (7) yields: 


n = ^- [ 
4 • jt J 


cos(a) 


2 


(9) 


This equation is in a form that is useful for solving for the downstream particle density when the emission site is the 
reservoir. When the emission site is a wall, it is more useful to re-express Equation (9) using a particle flux. This is 
done by solving Equation (6) for n w and substituting this into Equation (7) to obtain: 


n = JJjliU . f ( v ). d 3 v = J f(v ). dv . J ± 


4-r„ 


•dra n 


( 10 ) 


Note that T w appears with the second integral because the wall flux will be shown to be a function of the location on 
the wall. Substituting Equation (8) into Equation (10) yields: 


n = — • [ T w (A w ) 

C • 71 J 


cos(a) 


•dA„ 


(11) 


which is in a form that is useful for solving the particle density when a reflecting wall is the emission site. 


III. Equation for the Particle Flux From the Aperture Walls 

The particle fluxes and densities surrounding a cylindrical aperture are a function of particle fluxes from the 
reservoir and cylindrical aperture walls. While the particle flux from a reservoir is given by Equation (6), the wall 
flux has yet to be determined. The following section will review the development of the wall flux equation used in 
this study. 

Figure 2 is a schematic of the cylindrical aperture and the coordinate system used in this study. A cylindrical 
coordinate system will be used for the cylindrical aperture of radius R and thickness t. Also illustrated in Figure 1 
are the labels for the relevant surfaces and volumes. These labels will be used as subscripts within this study. It is 
convenient to non-dimensionalize the resulting equations, as was done in References 6 and 7. All dimensions will, 
therefore, be non-dimensionalized with the aperture radius, R, all particle fluxes will be non-dimensionalized with 
reservoir flux, T 0 , and all particle densities will be non-dimensionalized with reservoir density, n Q . Non- 
dimensionalized values will be denoted with a bar above the variable throughout this paper. 

Downstream 
Volume, “c” 



Upstream Volume 
or Reservoir, “o” 

Figure 2. — Schematic of the cylindrical aperture for this investigation. 
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The particle flux to the aperture walls has been solved many times before in literature (see, for example, Refs. 6, 
7, and 12). The following derivation can be found in References 6 and 7. The particle flux from the wall (i.e., 
surface “w”) of the aperture consists of the particle flux from upstream reservoir (i.e., surface “o”) arriving directly 
onto the wall and the particle flux from the other wall regions arriving back onto itself. Because of the symmetry of 
the problem, the particle flux from the wall is independent of tp. So, by applying Equation (5) to the reservoir and 
wall, the flux from the wall at any location z becomes: 


,W-— • | 

71 J 


cos(a 0 _ w )-cos(5 w _ 0 ) 


dA 0 + — • f T w (z) 

K J 


Ao 


cos(a w _ w )-cos(5 w _ w ) 


dA, 


( 12 ) 


Note here that a two digit subscript denotes two different surfaces. So, for example, L 0 _ w refers to the distance L 
between elemental areas on the wall and inlet plane (see Fig. 1). 

To solve this equation, cos(a), cos(5), L, and A for each integral must be solved as functions of z and the 
respective integral variable. The resulting solution in non-dimensional form is given by (Refs. 6 and 7): 


r w (z)=-- 


z 2 + 2 




z 2 + 4 


1 r w ( z w ) 


!-lz„-z|. 


(z w -z) 2 +4p 


' dz, 


(13) 


This equation has the form of a Fredholm integral equation of the second kind. This integral equation was solved 
numerically using an algorithm identical to that described in Reference 7. The resulting non-dimensional particle 
fluxes from the aperture wall are plotted as a function of the normalized wall axial location for various aperture 
thickness-to-diameter ratios in Figure 3. The figure illustrates a critical characteristic of free molecular flow through 
an aperture: the particle flux from an aperture wall as a function of wall axial location is linear over a broad range 
aperture thickness-to-diameter ratios. Indeed, as the aperture thickness-to-radius decreases, the functional relation 
becomes more linear. Over the range of thickness-to-radius ratios examined by this author in Figure 3, a linear 
assumption accurately models the wall particle flux as a function of wall axial location. Other characteristics 
include: 

1 . The wall particle flux as a function of wall axial location has a negative slope. This is because the source of 
particles is the upstream reservoir. 

2. The non-dimensional particle flux approaches 1 at the entrance and 0 at the exit as the thickness-to-radius 
ratio increases. Again, this is because the source of particles is the upstream reservoir. 

3. Over the range of thickness-to-radius ratios examined, the non-dimensional wall particle flux was 1/2 at the 
middle of the wall. 


Given the linear behavior of wall particle flux as a function of axial location, the integral equation of 
Equation (13) can be solved by assuming a linear equation for the solution. Clausing was the first to do this, and 
found the solution to be (Ref. 4): 


T w ( z ) = a + b • z 


where a and b are given by: 


+ 1 


4 


r +4 - 


t 2 


i 


r +4 +2 


and: 


respectively. 



(14) 


(15) 


(16) 
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Figure 3. — Non-dimensional particle flux from the aperture walls as a function of 
the normalized wall axial location for various aperture thickness-to-radius ratios. 
Aperture locations were normalized with the aperture thickness, t. 



Figure 4. — Comparison particle fluxes from the aperture walls as a function of 
the normalized wall axial location for various aperture thickness-to-radius 
ratios. Solid lines are numerically-determined and the dotted red lines are the 
respective values determined assuming a linear distribution. The wall axial 
location was normalized to the aperture thickness, t. 


Figure 4 compares the wall particle flux predicted by Equations (14) to (16) to that determined numerically. The 
worst-case difference between the two methods for aperture thickness-to-radius ratios up to 32 was less than 0.01 (in 
non-dimensional form). 

The remainder of this study will therefore assume that the wall particle flux as a function of wall axial location 
is given by Equations (14) to (16) because the assumption used to develop these equations is valid over the aperture 
thickness-to-radius ratios of interest to ion thruster grids (i.e., an aperture thickness-to-radius ratio of <1.25). The 
efficacy of this linear assumption will be assessed later in the results section. 
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IV. Equations for the Particle Flux and Densities 
Surrounding a Cylindrical Aperture 

In this section, the equations for determining particle fluxes and densities surrounding a cylindrical aperture will 
be determined. The first section develops the base equations. The second section divides the surrounding aperture in 
volumes, some of which are further divided into separate regions, and develops the limits of integration for these 
different cases. The last section presents the final equations for particle flux and density, as well as the transmission 
probability. 


A. Development of the Base Equations for Particle Flux and 
Density Surrounding a Cylindrical Aperture 

In this section, the base equations for determining the flux and density surrounding a cylindrical aperture are 
determined. Particle fluxes and densities surrounding an aperture are made up of contributions from the upstream 
reservoir and the aperture walls. Using Equation (5), the flux surrounding an aperture can be shown to be given by: 


r =—< 

n 


I 


cos(a 0 _ c ) • cos(5 0 _ c ) 


dA 0 + 


j r w (A w ) 


cos(a w _ c ) • cos(5 w _ c ) 


dA„ 


(17) 


Likewise, using Equations (9) and (11), the density surrounding an aperture can be shown to be given by: 



eos(a 0 -c) dA 


j r w (A w ) 


cos' 


(“w-c) 


• dA,, 


(18) 


Figure 5 shows a sketch of the geometric variables used in the equations above. For both equations, the left-hand 
integral is the reservoir’s contribution while the right-hand is the aperture wall’s contribution. Note that the reservoir 
flux and density are constants in Equations (17) and (18), respectively, while the wall contributions ofboth 
equations are a function of the wall location. And the origin of the coordinate system is the upstream center of the 
cylindrical aperture. Finally note that as shown in Figure 5, the flux determined in this study is in the axially 
downstream direction. 

To solve Equations (17) and (18), the variables a 0 _ c , 5 0 _ c , L 0 _ c , Ow_ c , 5 W _ C , and L w _ c and must be defined as 
functions of R, t, r 0 , r w , (p, r, and z. For the reservoir’s contribution, it can be shown that: 

L 0 -c = <Jz 2 +r2 +r o -2-r-r 0 -cos(<p) (19) 


and: 


a 0 _ c = 5 C _ 0 = arccos 


For the aperture wall’s contribution, it can be shown that: 

L,„_ 


■\j( z-z w ) 2 +r 2 +R 2 - 2-r-R-cos((p) 
R - r ■ cos((p)| 


and: 


5 C _ W = arccos 


(20) 


(21) 

(22) 


(23) 
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Downstream Volume 


Downstream Volume 



Upstream Reservoir 


Upstream Reservoir 


a. Emission from upstream reservoir. b. Emission from walls. 

Figure 5. — Sketch of the aperture geometry for determining particle fluxes and densities from an upstream reservoir 
and wall to a collection surface. Although the collection site is shown to be in the downstream region, the resulting 
equations will be applicable to all regions. 


Equations (19) to (23) can be substituted into Equations (17) and (18) to yield the equations for the flux and 
density surrounding a cylindrical aperture. Note that there is symmetry about <p, so numerical integration can be 
simplified by integrating only half of the cp integral and multiplying the result by 2. The resulting flux is given in 
non-dimensional form by: 


( J l, _V-_ 

[<p„j f 0 i [z “ + r + r 0 - 2 • r • r 0 • cos(<p )\ <p w . z w i [(: 


■' r w (z w ).(z-z w ).(l-r.cos(cp)) ^ I (24) 


J J /_ _ \2 -? i - 

(p wj z wi (z — z w j + r + 1 — 2 • r • c 


while the particle density is given by: 


J J 


TW,I “W, 

f o' d( P+ \ j 


r w( z w)-( 1 - r - cos ( ( p)) 


z +r +c - 2-r 


1— J 

• r 0 • cos(cp) 2 <p “ i Zwi (z-z w ) 2 + r 2 + 1 - 2 • r • cos(rp) 2 


-•dz w -d(pl (25) 


For both equations above, all dimensions (r, z, t, r 0 , and z w ) are non-dimensionalized with the aperture radius, R, all 
fluxes are non-dimensionalized with the reservoir flux, r o , and all densities are non-dimensionalized with the 
reservoir density, n 0 . The “i” and “f ’ subscripts for the limits of integration denote initial and final locations. It is 
possible to integrate with respect to r 0 and z w (i.e., the inner integral of each double integral) in the equations above 
after substituting Equation (14) into Equations (24) and (25). So, the equation for particle flux can be rewritten as: 

F(r,z)=r 0 (F,z)+r w (r,z) (26) 
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where the reservoir’s contribution is given by: 


r • r 0 • cos(cp)- z 2 - r 2 

r • cos(<p) 

— • urntnn 

r o - r • cos(<p) 


[r 2 ■ sin(<p)‘ + z 2 ]• [z 2 + r 2 + r 2 — 2 ■ r ■ r 0 

• C0S ( ( P)1 [f 2 -sin((p) 2 +z 2 f 

(XL U lull 

Jr 2 -sin((p) 2 +z 2 



—2 ^ 0 ’ f 

r 0 (i,z) = — • j 

K J 


and the wall’s contribution is given by: 

1 lPw,f \ 

r w (bz) = -- J [l-f-cos(cp)]- 


■ dtp 


(27) 


a + b-z„ 


(z-z w ) 2 +r 2 +l-2-r-cos(<p) J? 2 + 1 -2-r -cos(cp) 

Likewise, the equation for the particle density is given by: 

n(r, z )=n 0 (r, z)+n w (r, z) 
where the reservoir’s contribution is given by: 


<Po,f 


• arctan 


-m r"' f 


•Jr 2 + l-2-r-cos(cp) 


•d(p 


(28) 


n o( r ’ Z )“T“' [ 
2 • it J 


<Po,i 



r-r 0 -cos((p)-z 2 -r 


z 2 + r 2 • sin(cp ) 2 

'V 

z 2 + r 2 + rj — 2 • r • r 0 ■ cos((|)) 


•dtp 


and the wall’s contribution is given by: 


, <Pw,f 

— • f 

(r-cos((p)-l)- 

(a + b- z)- (z - z w ) + b • (r 2 + 1 - 2 • r • cos((p)) 

2 • n J 

<Pw,i 

[r 2 + 1 - 2 • r 

•cos((p) -^(z-z w ) 2 + r 2 + 1 - 2 • r • cos(tp) 


(29) 


(30) 


•dtp (31) 


B. Development of Limits of Integration for the Flux and Density Base Equations 

With the base equations for particle flux and density developed, all that remains is defining the limits of 
integration. However, the limits of integration are not necessarily constant but can be variable (i.e., functions of r, z, 
and the aperture geometry). As shown in the prior section, the particle flux and density at any location surrounding a 
cylindrical aperture is made up of particles that have a direct line-of-sight from the reservoir and the aperture walls. 
The limits of integration are, therefore, a function of this desired location (i.e., r and z) as well as the aperture 
geometry. Specifically, the contribution of the reservoir and walls depends on which portion of each area has a 
direct line-of-sight to the location of interest. This is illustrated in Figures 6 and 7. As shown in the figures, all 
reservoir and cylindrical wall areas have a direct line-of-sight to downstream locations within region I. But only a 
portion of the reservoir and cylindrical wall areas have a direct line-of-sight to downstream locations within 
regions II and III. And no reservoir particles have a direct-line of sight to downstream locations in region IV. 

In addition to this, there are three volumes surrounding a cylindrical aperture where the particle flux 
and density are desired. These volumes include the downstream, aperture interior, and upstream volumes. 
The equations developed for the downstream volume will be applicable to the interior volume, but a 
separate set of equations will be required for the upstream volume. 

Though complicated, the division of the space surrounding a cylindrical aperture into separate volumes and 
regions with different limits of integration illustrated in Figures 6 and 7 will used to develop a set of equations that 
can rigorously determine the particle flux and density surrounding a cylindrical aperture. The following sections will 
determine the borders between the different regions and their respective limits of integration. 
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Upstream 
Reservoir, n c 


Cylindrical 

Aperture 

Cross-section 


Downstream 

Volume 



Region I 



Shaded Reservoir 
Portions Have 
Direct Line-of-sight 
to Downstream 
Region 


Figure 6. — Definition of downstream regions for the development of the variable limits of 
integration for the particle flux from an upstream reservoir. 


Cylindrical 

Aperture 

Cross-section 


Upstream 
Reservoir, n 0 


1 


Region III 


> Region II 


> Region I 


Downstream 

Volume 


a. Downstream and Interior Volumes 





Shaded Wall 
Portions Have 
Direct Line-of-sight 
to Downstream 
Region 


CI> 




Region III 


Region II < 


Region I 


Shaded Wall 
Portions Have 
Direct Line-of-sight 
to Upstream Region 


Upstream 
Reservoir, n 0 


Cylindrical 

Aperture 

Cross-section 


A r 


b. Upstream Volume 


Downstream 

Volume 


Figure 7. — Definition of regions for the development of the variable limits of integration 
for the particle flux from the aperture walls. 
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1. Reservoir Contribution: Interior Volume and Downstream Volume — Region I 

As shown in Figure 6, there is a direct line of sight from any location to the entire reservoir surface as long as 
the following radial and axial border conditions are met: 


r < 1 and z >0 (32) 

As a result, the limits of integration for Equations (27) and (30) are constants (i.e., independent of r and z) and equal 
to 0 to jt for cp and 0 to 1 for r 0 . 

2. Reservoir Contribution: Downstream Volume — Region II 

Figure 6 shows region II of the downstream volume. It can easily be shown that region II is defined by the 
following radial and axial border conditions: 


1 < r < 4 and z > t (33) 

t 

Figure 8 shows the geometry of region II in greater detail. The shaded region in the figure has two different radii: r oII 
when cp is between 0 and <p 0 n; and R when <p is between (p oII and jt. As a result, integration for this case is slightly more 
complicated because the integrals of Equations (27) and (30) must be split into the summation of two integrals. For the 
first integral, it can be inferred from the figure that limits of integration along cp from 0 to <p oII will be between 0 and r oII 
for r 0 . The second integral’s limits of integration along (p between (p 0lI and n will be from 0 to R for r 0 . 

The limits r oII and (p oII must be defined as functions of r and z. The limit r oII is also a function of tp, which can be 
inferred from Figure 8, so the double integrals of Equations (27) and (30) must be integrated with respect to r 0 first. 
The limit r 0lI is determined by noting that triangles at axial locations t and z in Figure 8 share a common angle, y. 

The law of cosines is applied to the angle <p for the triangle at z to determine tan(a). The law of cosines is applied 
again to both triangles for the angle y and the cos(y) terms of each equation are equated to determine r oII . It can be 
shown that r on is given in non-dimensional for by: 



Upstream Reservoir 

Figure 8. — Development of the reservoir’s 
contribution to region II of the downstream 
volume. The shaded reservoir region 
contributes to the downstream flux and 
density. 
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To solve for <p 0 n, the equation above is set equal to R (or 1 in non-dimensional form) and solved for (p. The result in 
non-dimensional form is given by: 


1- 


r t \ 2 

I - 


(Poll (rt z ) = arccos 


r - _\2 
' t • r ' 

v z 


2* Vt .[ 1 - ‘ 


(35) 


3. Reservoir Contribution: Downstream Volume — Region III 

Figure 6 shows regions Ilia and Illb of the downstream volume. Region III is similar to region II, with the only 
difference being that shaded region in Figure 9 does not include the center of the aperture as it does in Figure 8 for 
region II. In fact, the equations for r 0lI and (p 0 n are still applicable in region III. However, the geometry of this region 
is more complicated, requiring that it be further divided into two separate regions (i.e., Ilia and Illb). The reason for 
this can be better understood by examining Figure 10, which shows the reservoir’s contribution at the aperture 
entrance. In this figure, the projected area is the circular area formed when viewing upstream from a downstream 
location through the aperture exit plane. It can be shown that the radius of this projected area, R p , is R/(l-t/z). The 
variable tp omax is the angle along tp where the radii of the aperture entrance and the projected area circles intersect to 
form a 90° angle, as shown in Figure 10. It can be shown that: 


. (r, z) = n - arcsin 

u - 


(36) 


Downstream Volume 



Upstream Reservoir 

Figure 9. — Development of the reservoir’s contribution to region III of the 
downstream volume. The shaded reservoir region contributes to the 
downstream flux and density. 
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Projected Area 


Projected Area 


a. Region Ilia b. Region Illb 

Figure 10. — Top view of the reservoir’s contribution to region III of the downstream volume. 

In region Ilia, the point of intersection for the radial lines used to define <p 0 max lies within the radius of the 
aperture, R. This creates some complexity for determining the limits of integration for Equations (27) and (30). For 
qK(p 0 n in region Ilia, it can be shown that the distance r 0 to the edge of the shaded region is r on as given in 
Equation (34). So, the limits of integration are from r 0lI to R between (p 0lI and n. However, between (p omax and cp 0 n> 
integration occurs between r oII and r oln , as shown in Figure 10. Using Figure 9 and the same process as that used to 
determine Equation (34), it can be shown that r oII i is given in non-dimensional form by: 


r oin( r >z,<p) = 




(37) 


In region Illb the point of intersection for the radial lines used to define (p 0 max lies outside the aperture entrance. 
This simplifies the the limits of integration for Equations (27) and (30). As shown in Figure 10, the distance r 0 to the 
edge of the shaded region is still r 0lI of Equation (34). And the limits of integration are from r oII to R between <p 0lI 
and 7i. 

The border between regions Ilia and Illb is determined by placing the point of intersection for the radial lines 
used to define <p 0 max at the aperture radius, R, and solving for the downstream radius as a function of axial location. 
So, region Ilia is defined by the following radial and axial border conditions: 


z 


— < r < 
t 


z 

t 



and z > t 


(38) 


Defining the border conditions for region Illb requires defining the border between regions Illb and IV. For 
Figure 6, it can easily be shown that this border is given by r = (2'z/t - R). So, region Illb is defined by the following 
radial and axial border conditions: 


z 


t 



+ 1 


<r< 



and z > t 


(39) 
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4. Reservoir Contribution: Downstream Volume — Regions IV 

As shown in Figure 6, there is no direct line of sight from the reservoir to region IV of the downstream volume. 
As a result, the particle flux and density contributions from the reservoir will be 0 for this region, which is defined 
by the following radial and axial border conditions: 


r > 



and z > t 


(40) 


5. Reservoir Contribution: Upstream Volume 

Narasimha solved for the particle flux and density surrounding a thin circular aperture (i.e., so thin that the 
thickness can be neglected) (Ref. 13). He found that there was symmetry about the aperture such that the upstream 
density could be determined with downstream density. For a spherical coordinate system, the reservoir particle 
density is given in non-dimensional form by: 


n o,u P ( r >7t-e) = l-n 0 (r,e) (41) 

where n oup is the reservoir’s contribution to the upstream particle density. Transforming coordinate systems from 
spherical to cylindrical, it can be shown that the reservoir’s contribution to the upstream particle density is given in 
non-dimensional form by: 


n o,up( r ,z) = l-n 0 (r,|z|) 


(42) 


Here, n 0 is given by Equation (30). By analogy, the reservoir’s contribution to the upstream particle flux is given in 
non-dimensional form by: 


r o, up ( r >z)=r o (r,|z|) 


(43) 


Here, T 0 is given by Equation (27). Later, it will be shown that Equations (27) and (30) decrease as z and r increase 
in magnitude. So at locations far upstream of the aperture, Equation (42) approaches n G (or unity in non-dimensional 
form), as would be expected because this is the reservoir’s density. However, Equation (43) approaches zero at these 
locations because the net axial flux should be zero at these locations due to the isotropic velocity assumption. 

6. Wall Contribution: Upstream, Interior, and Downstream Volumes — Region I 

As shown in Figure 7, there is a direct line of sight from any location to the entire aperture wall surface as long 
as the following radial border condition is met: 


r <1 (44) 

This is true for any volume (upstream, interior, or downstream), so there is no axial border condition. As a result, the 
limits of integration for Equations (28) and (31) are constants (i.e., independent of r and z) and equal to 0 to n for (p 
and 0 to t for z w . 

7. Wall Contribution: Downstream Volume — Region II 

Figure 7 shows region II of the downstream volume. It can easily be shown that region II is defined by the 
following radial and axial border conditions: 


1 < r < 



and z > t 


(45) 
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Downstream Volume 



Upstream Reservoir 

Figure 1 1 . — Development of the aperture wall’s contribution to 
region II of the downstream volume. Only the shaded wall 
surface contributes to the downstream flux and density. Only 
half of the shaded region is shown. 


Figure 1 1 shows the geometry of region II in greater detail. Integration for this case is slightly more 
complicated in that the integrals of Equations (28) and (31) have to be split into the summation of two integrals. It 
can be inferred from the figure that the first integral’s limits of integration will be between (p wt and <p w n, for z w from 
z wII to t. The second integral’s limits of integration will be between <p wII and n for z w from 0 to t. 

The limits z wII , <p w n, and <Pwt must be defined as functions of r and z. The limit z wII is also a function of cp, which 
can be inferred from Figure 11, so the double integrals of Equations (28) and (31) must be integrated with respect to 
z w first. The limit z wII is determined as was done for the reservoir. Note that the triangles at z = t and z in Figure 1 1 
share a common angle, y. The law of cosines is applied to the angle (p for the triangle at z to determine tan(P). The 
law of cosines is applied again to both triangles at t and z for the angle y and the cos(y) terms of each equation are 
equated to determine z wII . It can be shown that z wn is given in non-dimensional for by: 


/ r 2 ■ t + t - 2 • z - 2 • r • (t - z)- cos(< 

z wli V r > z > 9 1 = 


F 2 -! 


(46) 


To solve for <p wII , the equation above is set equal to 0 and solved for cp. The result in non-dimensional form is given 
by: 


(p wII (r,z) = arccos 


r 2 +l-24 


2 -f • 1- 


(47) 


To solve for (p wt , the equation for z wII above is set equal to t and solved for (p. The result in non-dimensional form is 
given by: 


cp wt (r,z)=arccos| ^ 


(48) 
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8. Wall Contribution: Downstream Volume — Region III 

Figure 7 shows region III of the downstream volume. It can easily be shown that region III is defined by the 
following radial and axial border conditions: 


r > ^2 • 4 - 1 j and z > t (49) 

Figure 12 shows the geometry of region III in greater detail. Integration of Equations (28) and (31) for this case is less 
complicated than region II. It can be inferred from the figure that the integral’s limits of integration will be between <p wl 
and 7i for z w from z wD to t. Note in Figure 12 that the equations for cp*, and z w[I are the same as those for region II. 

9. Wall Contribution: Upstream Volume — Region II 

Figure 7 shows region II of the upstream reservoir. It can easily be shown that region II is defined by the 
following radial and axial border conditions: 


1 < r < 



and z < 0 


(50) 


Figure 13 shows the geometry of region II in greater detail. As with the downstream case, integration for this case is 
slightly more complicated because the integrals of Equations (28) and (31) have to be split into the summation of 
two integrals. It can be inferred from the figure that the first integral’s limits of integration will be between <p wt and 
(Pwiim where integration along z w will be between 0 and z wIIn . The second integral’s limits of integration will be 
between <p wIIn and 7t, where integration along z w will be between 0 and t. 


Downstream Volume 



Upstream Reservoir 

Figure 12. — Development of the aperture wall’s contribution to 
region III of the downstream volume. Only the shaded wall 
surface contributes to the downstream flux and density. Only 
half of the shaded region is shown. 
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Downstream Volume 



Figure 13. — Development of the aperture wall’s contribution to region II of the 
upstream volume. Only the shaded wall surface contributes to the downstream 
flux and density. Only half of the shaded region is shown. 


The limits z wIIn , <p w iin, and <Pwt must be defined as functions of r and z. The limit z wIIn is also a function of <p, 
which can be inferred from Figure 13, so the double integrals of Equations (28) and (31) must be integrated with 
respect to z w first. The limit z wIIn is determined as was done for the reservoir. Note that the triangles at z = 0 and -z in 
Figure 13 share a common angle, y. The law of cosines is applied to the angle (p for the triangle at z to determine 
tan(P). The law of cosines is applied again to both triangles at t and z for the angle y and the cos(y) terms of each 
equation are equated to determine z wIIn - It can be shown that z wIIn is given in non-dimensional for by: 


Zwiin (r, z,cp) = 2-z- 


l-r-cos((p) 


1-r 2 


(51) 


To solve for (p w n n , the equation above is set equal to t and solved for (p. The result in non-dimensional form is given 
by: 


<Pwlln(r,z) 


arccos 




(52) 


The equation for cp wt is given in Equation (48). 

10. Wall Contribution: Upstream Volume — Region III 

Figure 7 shows region III of the upstream volume. It can easily be shown that region III is defined by the 
following radial and axial border conditions: 


r > 



and z < 0 


(53) 


Figure 14 shows the geometry of region III in greater detail. Integration of Equations (28) and (31) for this case is 
less complicated than region II. It can be inferred from the figure that the integral’s limits of integration will be 
between cp wt and n, where integration along z w will be between 0 and z wn . Note in Figure 14 that the equations for 
cp wt and z wIIn are the same as those for region II. 
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Downstream Volume 



Figure 14. — Development of the aperture wall’s contribution 
to region III of the upstream volume. Only the shaded wall 
surface contributes to the downstream flux and density. 

Only half of the shaded region is shown. 

C. Equations for Flux and Particle Density Surrounding a Cylindrical Aperture 

The limits of integration defined in the previous section can be substituted into Equations (27) and (28) for 
particle flux and Equations (30) and (31) for particle density. Because the surrounding aperture had to be divided 
into separate volumes and regions as illustrated in Figures 6 and 7, the radial and axial border conditions defined in 
the previous sections are used to determine the appropriate particle flux or density equation. The results are listed in 
Tables 1 to 5. The particle flux in the axial direction is given by Equation (26), with the reservoir’s contribution 
equations listed in Table 1 and the wall’s contribution equations listed in Table 2. The particle density is given by 
Equation (29), with the reservoir’s contribution equations listed in Table 3 and the wall’s contribution equations 
listed in Table 4. Table 5 lists the variable limits of integration equations used in Tables 1 to 4. 

The integrals remaining in Equations (27), (28), (30), and (31) could only be solved in two cases, which 
included the flux equations for region I of the downstream and interior volumes and for the upstream volume 
(see Table 1). The remaining integrals in Tables 1 to 4 must be solved numerically. 

It can be shown that the particle density equation for the reservoir’s contribution in region I produces a 
singularity at r = z = 0. This is easily overcome by noting that the reservoir’s contribution to particle density is 1/2 at 
the aperture entrance (i.e., z = 0), and is independent of radial location (Ref. 11). 

The particle transmission probability, P, through the cylindrical aperture (i.e., the fraction of particles incident 
upon the aperture inlet that eventually pass through) can also be determined (Ref. 7). The probability is given by 
integrating the total particle flux radially at any axial location and dividing by the particle flux at the aperture inlet: 

2-71 GO 

| | r(r, z) • r • dr • dcp ^ 

p (z)= °° 2 -nR = 2-jF(r,z)-r-dr ( 54 ) 

r o ' j J r • dr • d <P ° 

o o 

Within the cylindrical aperture (i.e., 0 < z < t), the upper limit of integration along r changes from oo to R (or 1 in 
non-dimensional form). Because the transmission probability can be calculated at any axial location and still yield 
the same result, this equation also serves as a useful tool to verify the equations in Tables 1, 2, and 5. And because 
particle density equations utilize the same equations except for the cos(5) terms, it can indirectly verify the equations 
in Tables 3 and 4. 
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Upstream volume 





• arctan 




TABLE 3.— EQUATIONS FOR THE RESERVOIR’S CONTRIBUTION TO THE PARTICLE DENSITY FOR THE DIFFERENT REGIONS OF FIGURE 6 

Equation | Border conditionals 

Downstream and Interior volumes 



Upstream volume 







TABLE 5.— VARIABLE LIMITS OF INTEGRATION USED IN TABLES 1 TO 4 
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V. Far-fleld Equations Using Clausing’s Angular Distribution Equation 

A commonly accepted far-field expression will be presented here for later comparisons with the more rigorously 
developed equations above. Clausing was the first to develop a far-field expression for the angular distribution of 
particles downstream of a cylindrical aperture (Ref. 3). His equation assumes that the aperture is a point source and that 
the flux of particles from the walls of the aperture is linear as a function of axial thickness. The angular distribution of 
particles downstream of a cylindrical aperture, T(a), in the far-field is given by (Refs. 4 and 5): 


T(a) = 




„i n | 1 ' tan ( a ) j , I ' tan ( a ) [ 1 1 ' tan ( a ) 


1- 


" + ' (l _ 2 • k) • 

3 • 7T 


j [' ttan(a) 
2 


t • tan(a) 


( 2 

if a < arctan — 

l t 


(55) 


and 


4 1-2-k . 


K H 7 — r- 

3 * 7i t • tan (a) 


if a > arctan 


li 


where: 


4 


t 2 +4 -t 


(56) 


2 + 




/ r +4 

Figure 15 shows the geometiy for a single aperture that is treated as a point source. It can easily be shown that: 


L = 


■Jr 2 +(z-t) 2 


(57) 


and: 


cos(a) = cos(§)= Z { = (58) 

V'-Mz-O 2 

Equations (57) and (58), along with T(a), can be substituted into Equations (5) and (9) to determine the downstream 
particle flux and density from a cylindrical aperture in the far-field. Note that all variables are independent of dA 0 
because a point source is assumed, so that integral reduces to jtR 2 . The resulting equation for the far- field flux, T far , 
is given in non-dimensional form by: 


Tfar(r, z )= T 


arccos 


z- 1 


while the far-field density, nf ar , is given in non-dimensional by: 


M ) 2 

(f 2 +(z-tf 


n far(rt z ) = ^-T 


r ^ 

z - 1 


arccos 


A /r 2 +(z-t) 2 


z- 1 


2 ' (z-t ) 2 


n lr + 


(59) 


(60) 
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Figure 15. — Geometry for determining downstream 
axial flux and density in the far-field. 


The transmission probability of particles through a cylindrical aperture can also be determined by integrating 
Clausing’s angular distribution function over all azimuthal and polar angles and dividing by that of a diffuse 
distribution. It can be shown that this probability, P far , is given by (Ref. 5): 


TC 

2 -n 2 

J jT(a)-cos(a)-sin(a)-da-d(p ^ 

P far =-^-5 — = 2- jT(a)-cos(a)-sin(a)-da (61) 

2-jt 2 0 

J Jcos(a)-sin(a)-da-d(p 
0 0 


VI. Results 

The following sections will present the results of the particle flux and density equations developed in this study. 
First, an assessment of the linear wall flux distribution assumption will be made. Then, the axial flux and density 
equations from this study will be applied to a cylindrical aperture and compared to those of Clausing’s far-field 
expression to accomplish the following: 1) the commonly accepted far-field expression will be used to verify the 
rigorously developed equations in the far-field; and 2) the rigorously developed equations will be used to determine 
the efficacy of the far-field expression in the near-field. A typical ion thruster aperture thickness-to-radius ratio of 
1.25 will be used throughout. Finally, the axial flux and density equations will be verified using three different 
methods. 
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A. Assessment of Linear Wall Flux Assumption and Transmission Probabilities 

As was shown earlier in Figures 3 and 4, the linear axial wall flux distribution assumption of Equation (14) 
accurately models the true distribution over a range of aperture thickness-to-radius ratios. However, it would be 
useful to determine a figure of merit for setting the maximum thickness-to-radius ratio that can utilize this 
assumption. One figure of merit is the particle transmission probability. Figure 16 shows the error between the true 
transmission probabilities (i.e., those determined numerically per Reference 7 with an error of <0.25 percent) and 
those determined with Equation (54), which assumes a linear distribution assumption. The figure shows that 
transmission probability errors with a linear wall flux distribution assumption are 4 percent for thickness-to-radius 
ratios <16. However, in state-of-the-art ion thruster perforated grids, the thickness-to-radius ratio is within 1.25, so 
the transmission probability error would be <0.1 percent. 

Verifying the linear distribution assumption as the cause for this error is essential for assessing its impact on axial 
flux and density calculations. To better understand the cause of the transmission probability errors, Figure 17 plots the 
transmission probability of Equation (54) as a function of axial location for a thickness-to-radius ratio of 32. This 
thickness-to-radius ratio was selected because these errors increase at larger ratios. The transmission probabilities 
should be constant throughout the axial direction in Figure 17. However, the figure shows that the probability is not 
constant within the aperture. So, the transmission probability can be in error by as much as 1 1 percent for this 
thickness-to-radius ratio. It can be shown that this disparity generally increases for larger thickness-to-radius ratios. 
Also shown in Figure 17 is the transmission probability when the linear wall flux distribution assumption is replaced 
with a 5 th order polynomial curve fit equation. The polynomial curve fit equation produces a nearly constant 
transmission probability as a function of axial location whose average value is nearly identical to the true transmission 
probability. This is significant, because the flux difference between the hue and linear distributions never exceeded 
1 percent ofT 0 , and yet the impact on transmission probability was a nearly 7 percent error. 



0 4 8 12 16 20 24 28 32 

Thickness-to-Radius Ratio, F 

Figure 16. — True transmission probability as a function of aperture thickness-to- 
radius ratio and the error between the true probability and that assuming a 
linear wall flux distribution. 


NASA/TM— 20 11-216970 


27 



0.076 


0.074 


0.072 


0.07 


0.068 



Non-dimensional Axial Locations 


Figure 17. — Transmission probability with linear and polynomial curve fits to the 
actual wall flux distribution for a thickness-to-radius ratio of 32. 


Because the transmission probability is the integral of the axial flux per Equation (54), the errors plotted in 
Figure 16 would roughly represent the errors expected for flux, and therefore density, calculations. This assertion 
was checked for a thickness-to-radius ratio of 32. While the maximum transmission probability error was about 
7 percent, the maximum flux error was within 8 percent. So, if a maximum acceptable flux and density error is 
4 percent, a linear wall flux distribution assumption can be used for thickness-to-radius ratios <16 according to 
Figure 16. 


B, Particle Flux and Density Surrounding a Cylindrical Aperture 

1. Near- field Results 

The axial flux surrounding a cylindrical aperture with a thickness-to-radius ratio of 1.25 is shown in Figure 18. 
The plot was determined with Equation (26) and the equations from Tables 1, 2, and 5. As expected, the axial flux 
decreases rapidly further from the aperture in the upstream, downstream, and radial directions. This is because the 
flux tends to be equal in all directions upstream of the aperture and the flow expands rapidly downstream of the 
aperture. It is only within the cylindrical aperture that the aperture and its walls align the flux of the particle flow in 
the axial direction. It is interesting to note that there is symmetry in the axial direction about z = t/2. 

Figure 19 shows the contributions of the reservoir and aperture walls to the surrounding axial flux. As expected, 
the non-dimensional axial flux from the reservoir is 1 (or the reservoir flux, r o ) at the aperture inlet plane and this 
decreases rapidly as the flow expands in all downstream directions. The wall’s contribution produces a net negative 
(i.e., upstream) axial flux near the aperture inlet because of upstream particle reflections from the aperture walls. 

The particle density surrounding a cylindrical aperture with a thickness-to-radius ratio of 1 .25 is shown in 
Figure 20. The plot was determined with Equation (29) and the equations from Tables 3 to 5. Far upstream from the 
aperture inlet, the non-dimensional particle density is constant and equal to 1 (or the reservoir density, n 0 ). Flowever, 
closer to the aperture inlet, the figure shows a depleted density zone. This depletion zone is caused by particles 
escaping through the aperture. The figure also shows that the density drops rapidly from the upstream to 
downstream regions and through the aperture. Figure 21 shows only the particle density downstream of the aperture. 
The particle density drops rapidly as particles expand into vacuum, dropping to less than 10 percent of the upstream 
density within 1 aperture radius of the downstream surface. 
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Figure 18. — Axial particle flux surrounding a cylindrical aperture 
with a thickness-to-radius ratio of 1.25. 
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a. Reservoir’s contribution. b. Wall’s contribution. 

Figure 19. — The reservoir’s and wall’s contributions to the axial particle flux surrounding a cylindrical aperture 
with a thickness-to-radius ratio of 1 .25. 
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Figure 20. — Particle density surrounding a cylindrical aperture 
with a thickness-to-radius ratio of 1 .25. 



Figure 21 . — Particle density downstream of a cylindrical 
aperture exit plane with a thickness-to-radius ratio of 
1.25. 
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2. Far-field Results and Comparisons with Clausing ’s Far-field Expression 

The far-field particle flux and density equations (i.e., Eqs. (59) and (60), respectively) that utilize Clausing’s 
angular distribution function are compared to that calculated with Equations (26) and (29), respectively, and the 
equations of Tables 1 to 5 in Figures 22 and 23. Here, the axial flux and density were determined at fixed radial 
distances from the downstream center of the aperture from 2 to 60 radii while the angle was varied (i.e., a spherical 
coordinate system was used) to determine angular distributions. The figures illustrate that at increasing distances from 
the aperture, the far-field angular distribution becomes equal to that of the Equations (26) and (29). The figures also 
show that the worst-case error in the far-field equation occurs along the axial centerline of the aperture. 
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Figure 22. — Normalized axial flux angular distributions for a variety of radial 
distances in a spherical coordinate system. The aperture thickness-to- 
radius is 1.25. Solid lines were determined from Equation (59) and dashed 
lines from Equation (26) and Tables 1, 2, and 5. 
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Figure 23. — Normalized density angular distributions for a variety of radial 
distances in a spherical coordinate system. The aperture thickness-to- 
radius is 1 .25. Solid lines were determined from Equation (60) and 
dashed lines from Equation (29) and Tables 3 to 5. 
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Because Equations (26) and (29) along with the equations of Tables 1 to 5 represent a rigorous solution for 
downstream particle flux and density, respectively, the far-field flux and density equations were compared to them 
in Figures 24 and 25 to determine the error as a function of radial distance from the downstream aperture center at 
various angles. The figures show that both far-field equation errors reduce to about 5 percent at about 25 aperture 
radii downstream for the 0° case, but within 9 aperture radii for the remaining angles. These figures demonstrate that 
the far-field equations will produce accurate results at >9 aperture radii downstream of a cylindrical aperture with a 
thickness-to-radius ratio of 1.25. 



Non-dimensionalized Radial Distance 

Figure 24. — Far-field flux error as a function of radial distance at various 
angles in a spherical coordinate system. The aperture thickness-to-radius 
is 1.25. Error is equal to 100 |Equations (26) to (59)|/Equation (26). 
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Figure 25. — Far-field density error as a function of radial distance at various 
angles in a spherical coordinate system. The aperture thickness-to-radius is 
1.25. Error is equal to 100 |Equations (29) to (60)|/Equation (29). 


NASA/TM- 


-2011-216970 


32 



C. Verification of Equations — Clausing’s Far-field Expression and Conservation of Mass 

The equations determined in this study for the axial flux and density surrounding a cylindrical aperture were 
verified using three methods. First, comparisons were made to a commonly accepted expression determined by 
Clausing (Refs. 3 and 4). This was done in the prior section for both the axial flux and density in Figures 22 to 25. 
Results showed that maximum differences were within 5 percent at >25 aperture radii downstream and that these 
differences decreased with increasing distance. So, the equations determined in this study compare favorably to that 
determined by Clausing in the far-field downstream volume. 

A second method was to compare the transmission probabilities of the far-field expression (i.e., Eq. (61)) to that 
of the more rigorous Equation (54). Because both equations were determined assuming a linear wall flux 
distribution, both equations should yield the same transmission probabilities. Comparisons were made over a broad 
range of thickness-to-radius ratios, and Equations (54) and (61) were found produce nearly identical results. 

A third method was to verify conservation of mass throughout the axial direction. For mass to be conserved, the 
total axial particle flow rate must be constant (i.e., independent of axial location, z). The total axial particle flow rate 
is determined by integrating the axial flux radially at a given location. Note that this is merely the transmission 
probability multiplied by the particle flux at the aperture inlet, which is a constant. Therefore, merely demonstrating 
a constant transmission probability along z with Equation (54) would assure that mass was conserved. This was 
done over a range of thickness-to-radius ratios and axial locations, and transmission probabilities were found to be 
constant, so mass was conserved. 


VII. Conclusions 

The equations for rigorously calculating the particle flux and density surrounding a cylindrical aperture in the 
free molecular flow regime were developed and presented. Previous work did not provide accurate near-field 
determinations of particle flux and density. The fundamental equations for particle flux and density from a reservoir 
and a diffusely reflecting wall were developed and presented. Assumptions included, in part, a Maxwell-Boltzmann 
speed distribution, equal particle and wall temperatures, and a linear flux distribution from the cylindrical aperture 
walls. With this information, the equations for axial flux and density surrounding a cylindrical aperture were 
developed. The cylindrical aperture was divided into three volumes (i.e., upstream reservoir, interior, and 
downstream), and the upstream and downstream volumes were further divided into multiple regions to rigorously 
determine axial flux and density. After the base equations for flux and density were determined, the limits of 
integration were developed for all volumes and regions. The complete set of flux and density equations was 
presented along with an equation for the transmission probability. Finally, a commonly accepted far-field angular 
distribution expression was used to develop far-field flux, density, and transmission probabilities for comparisons. 

An assessment of the linear wall flux distribution assumption showed that a thickness-to-radius ratio of 16 
would produce an axial flux and density error of 4 percent and that this error decreased with decreasing thickness-to- 
radius ratios. The axial flux and density surrounding a cylindrical aperture with a thickness-to-radius ratio of 1 .25 
was presented. The axial flux was found to be axially symmetric about 1/2 of the aperture plate thickness. Both axial 
flux and density were found to decrease rapidly downstream as the flow exhausting from the aperture expanded into 
vacuum. 

Finally, the equations determined in this study were verified using three methods. Comparisons to a commonly 
accepted far-field expression showed that maximum differences were within 5 percent at 25 aperture radii 
downstream and that these differences decreased with increasing distance. The transmission probabilities of the far- 
field expression were also nearly identical to those determined from this study. Finally, conservation of mass was 
confirmed in the axial direction. 
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